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Abstract. Vacuum correlator of topological charge densities x(<Z 2 ) a t low q 2 
£T) \ in QCD is discussed: its value x(0) and first derivative x'(0) at q 2 = and the 

C\l ■ contributions of pseudoscalar quasi- Goldstone bosons to x(l 2 ) a t l° w Q 2 - The 

(N ; qcd sum rule 111] , giving the connection of x'(0) (f° r massless quarks) with 

the part of the proton spin E, carried by u, d, s quarks, is presented. From the 
requirement of selfconsistancy of the sum rule the values of x'(0) and S were 
■ found. The same value of x'(0) follows also from the experimental data on 

S. The contributions of 7r and r\ to x(q 2 ) are calculated basing on low energy 
theorems. In such calculation the ir — rj mixing, expressed in terms of quark 
mass ratios is of importance. 
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1. Introduction. 



■ The existence of topological quantum number is a very specific feature of non- 



abelian quantum field theories and, particularly, QCD. Therefore, the study of 
properties of the topological charge density operator in QCD 



8tt 

and of the corresponding vacuum correlator 



sO*) = ^Gl v {x)Gl v (x) (1) 



X (q 2 ) = i / d A xe^(i) | T {Q 5 (x), Q 5 (0)} | 0) (2) 



is of a great theoretical interest. (Here G™„ is gluonic field strength tensor, G^ v = 
(l/2)e ^xaGxa is its dual, n are the colour indeces, n — 1,2, ...N 2 — 1, N c is the 
number of colours, _/V c = 3 in QCD). The existence of topological quantum numbers 
in non-abelian field theories was first discovered by Belavin et al. , their connection 
with non-conservation of U(l) chirality was established by t'Hooft p|. Gribov B 
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was the first, who understood, that instanton configurations in Minkowski space 
realize the tunneling transitions between states with different topological numbers. 
Crewther |5| derived Ward identities related to x(0), which allowed him to prove 
the theorem, that x(0) = in any theory where it is at least one massless quark. 
An important step in the investigation of the properties of x(<7 2 ) was achieved by 
Veneziano |(| and Di Vecchia and Veneziano M . These authors considered the limit 
N c — ► oo. Assuming that in the theory there are Nf light quarks with the masses 
in i <C M, where M is the characteristic scale of strong interaction, Di Vecchia and 
Veneziano found that 



X(0) = <0|gg|0)^ ±) , (3) 

where (0 | qq | 0) is the common value of quark condensate for all light quarks 
and the terms of the order m,/M are neglected. [] The concept of 6>-term in 
the Lagrangian was succesfully exploited in [g] in deriving of (|^). Using the same 
concept and studying the properties of the Dirac operator Leutwyler and Smilga || 
succeeded in proving eq.3 at any N c for the case of two light quarks, u and d. How 
this method can be generalized to the case of Nf = 3 is explained in the recent 
Paper §. 

In this paper I discuss x(<? 2 ) in the domain of low | q 2 |, i.e. I suppose that 
q 2 /M 2 is a small parameter and restrict myself to the terms linear in this ratio. In 
this domain must be accounted terms x(0), x'(0)<Z 2 as well as the contributions of 
low mass pseudoscalar quasi-Goldstone bosons it and T), resulting to nonlinear q 2 
dependence. 

Strictly speaking, only nonperturbative part of x(g 2 ) has definite meaning. The 
perturbative part is divergent and its contribution depends on the renormalization 
procedure. For this reason only nonperturbative part of x(<? 2 ) @ with perturbative 
part subtracted will be considered here. In ref.l on the basis of QCD sum rules in 
the external fields the connection of x'(0) (its nonperturbative part) with the part of 
the proton spin E, carried by u, d, s quarks, was established. From the requirement 
of the selfconsistancy of the sum rule it was obtained 

X'(0) = (2.4 ± 0.6) x 1(T 3 GeV 2 (4) 

in the limit of massless u, d and s quarks. The close to (f|) value follows also 
by the use of experimental data on E. At low q 2 the terms, proportional to quark 
masses, are related to the contributions of light pseudoscalar mesons as intermediate 
states in the the correlator (||). These contributions are calculated below. In such 
calculation for the case of three quarks the mixing of tt° and 77 is of importance and 
it is accounted. 

The presentation of the material in the paper is the following. In Sec. 2 the QCD 
sum rule for E is derived. The important contribution to the sum rule comes from 



1 The definition of x(l 2 ) used above in eq.(|^), differs by sign from the definition used in Q-^]- 
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x'(0). This allows to determine x'(0) m two ways: using the experimental data on £ 
and from the requirement of the selfconsistancy of the sum rule. Both methods give 
the same value of x'(0)> given by eq.[|. In Sec. 3 the low energy theorems related to 
x(0) arc rederived with the account of possible anomalous equal-time commutator 
terms. (In ||-[[7| it was implicitly assumed that these terms are zero). In Sec. 4 the 
case of one and two light quarks are considered. It is proved, that the mentioned 
above commutator terms are zero indeed and for the case of two quarks eq.3 is 
reproduced without using N c — > oo limit and the concepts of 0-terms. In Sec. 5 the 
case of three it, d, s light quarks is considered in the approximation m u , m-d m s . 
The problem of mixing of 7r° and r/ states JTo|, |TT[] is formulated and corresponding 
formulae are presented. The account of tt — rj mixing allows one to get eq. (||) from 
low energy theorems, formulated in Sec. 2 for the case of three light quarks. (At 
m u , rrid <C m s it coincides with the two quark case). In Sec. 5 the g 2 -dependence of 
x(<? 2 ) was found at low | q 2 | in the leading nonvanishing order in q 2 /M 2 as well as 
in m q /M. 



2. QCD sum rule for S. Connection of S and x'(O). 

As well known, the parts of the nucleon spin carried by it, d and s-quarks are 
determined from the measurements of the first moment of spin dependent nucleon 
structure function gi(x,Q 2 ) 

l 

r p ,„(Q 2 ) = J dxg\- Ptn (x, Q 2 ) (5) 
o 

The data allows one to find the value of S - the part of nucleon spin carried by 
three flavours of light quarks S = Au + Ad + As, where Am, Ad, As are the parts 
of nucleon spin carried by u, d, s quarks. On the basis of the operator product 
expansion (OPE) S is related to the proton matrix element of the flavour singlet 
axial current j® 5 

2ms yY, = (p,s\j° 5 \p,s), (6) 

where s^ is the proton spin 4-vector, m is the proton mass. The renormalization 
scheme in the calculation of perturbative QCD corrections to T p . n can be arranged 
in such a way that S is scale independent. 

An attempt to calculate X using QCD sum rules in external fields was done in 
ref . |Q . Let us shortly recall the idea. The polarization operator 



n(p)=i / d 4 xe lpx (0\T{r]{x),fi(0)}\0) 



(7) 
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was considered, where 

r)(x) = e abc ^(xp^ix)] W d c {x) (8) 

is the current with proton quantum numbers |l3} |, u a , d b are quark helds, a,b,c are 
colour indeces. It is assumed that the term 

AL = 3%A, (9) 

where is a constant singlet axial field, is added to QCD Lagrangian. In the weak 
axial field approximation H(p) has the form 

n(p) = n(°)(p) + n«(p)A M . (io) 

is calculated in QCD by OPE at p 2 < 0, |p 2 | > -R^ 2 , where i? c is the confine- 
ment radius. On the other hand, using dispersion relation, Ilj} (p) is represented by 
the contribution of the physical states, the lowest of which is the proton state. The 
contribution of excited states is approximated as a continuum and suppressed by 
the Borel transformation. The desired answer is obtained by equalling of these two 
representations. This procedure can be applied to any Lorenz structure of (p) , 
but as was argued in [[l4j , [jl5| the best accuracy can be obtained by considering the 
chirality conserving structure 2p fl p"f5. 

An essential ingredient of the method is the appearance of induced by the 
external field vacuum expectation values (v.e.v). The most important of them in 
the problem at hand is 

(0\j%\0) A =3f^ (11) 

of dimension 3. The constant /q is related to QCD topological susceptibility. Using 
(^), we can write 

(0\j%\0)a = Um^o i J d 4 xe^(0\T{j° 5 (x),j%(0)}\0)A„ = 

= lim q ^ P^ u {q)A u (12) 
The general structure of P^iq) is 

P^{q) = -P L (q 2 )S^ + P T (q 2 )(-S^q 2 + q„q v ) (13) 

Because of anomaly there are no massless states in the spectrum of the singlet po- 
larization operator P^ even for massless quarks. Pr,i(9 2 ) also have no kinematical 
singularities at q 2 = . Therefore, the nonvanishing value P^{0) comes entirely 
from Pi(q 2 ). Multiplying P t _ LU {q) by q^q^, in the limit of massless u, d, s quarks we 
get 
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q»q»P»u{q) = -P L (q 2 )q 2 = N 2 (a s /4n) 2 i I d 4 xe l ^x 

x (0\TG^(x)G^(x),GZ(0)GZ(0)\0), (14) 

where G" is the gluonic field strength, G A ,„ — (l/2)e AI ^AcrGA ( T-(The anomaly con- 
dition was used, Nf = 3.). Going to the limit q 2 — > 0, we have 

/ 2 = -(1/3)P £ (0) = ^Y(O), (15) 

where x(<Z 2 ) is defined by (||). 

According to Crewther theorem ||), x(0) = if there is at least one massless 
quark. The attempt to find x'(0) itself by QCD sum rules failed: it was found Jl^j 
that OPE does not converge in the domain of characteristic scales for this problem. 
However, it was possible to derive the sum rule, expressing £ in terms of f 2 ( |Ti"| ) 
or x'(0). The OPE up to dimension d — 7 was performed in rcf. |l2|. Among the 
induced by the external field v.e.v.'s besides ([ll]), the v.e.v. of the dimension 5 
operator 

g(0\Y / ^(l/2)X n G^q\0)A = 3h Q Ap, q = u,d,s (16) 

was accounted and the constant ho was estimated using a special sum rule, 
ho ~ 3 x 10~ 4 GeV 4 . There were also accounted the gluonic condensate d = 4 and 
the square of quark condensate d = 6 (both times the external field operator, 
d = 1). However, the accuracy of the calculation was not good enough for reliable 
calculation of £ in terms of /q : the necessary requirement of the method - the weak 
dependence of the result on the Borel parameter was not well satisfied. 

In § the accuracy of the calculation was improved by going to higher order 
terms in OPE up to dimension 9 operators. Under the assumption of factorization 
- the saturation of the product of four-quark operators by the contribution of an 
intermediate vacuum state - the dimension 8 v.e.v.'s are accounted (times A^): 

- g(0\qa aP (l/2)X n G^q ■ qq\0) = m 2 (0\qq\0) 2 , (17) 

where m 2 = 0.8 ± 0.2 GeV 2 was determined in |l6|. In the framework of the same 
factorization hypothesis the induced by the external field v.e.v. of dimension 9 

^(OI^IOU^IO) 2 (18) 

is also accounted. In the calculation the following expression for the quark Green 
function in the constant external axial field was used |^|| : 

(0\T{q a a (x), $(0)}|0) A = i5 ah x aP /2K 2 x 4 + 



+ {l/2n 2 )8 ab {Ax)( l5 x) aP /x 4 - {l/\2)5 ah 5 a0 (Q\qq\O) + 



G 



-(l/72^ afc (0|gg|0)(ii 75 - Ax lb ) af3 - 



+ (1/ 12) f 2 8 ab (A l5 ) a p + (l/216)rt 



(5/2)x 2 ^7 5 - (Ax)x-y 5 



(19) 

The terms of the third power in x-expansion of quark propagator proportional to 
Ay, are omitted in (|l9]), because they do not contribute to the tensor structure of 
Tly of interest. Quarks are considered to be in the constant external gluonic field 
and quark and gluon QCD equations of motion are exploited (the related formulae 
are given in jL7|). There is also an another source of v.e.v. ho to appear besides 
the x-expansion of quark propagator given in eq.([l9|): the quarks in the condensate 
absorb the soft gluonic field emitted by other quark. A similar situation takes place 
also in the calculation of the v.e.v. ( |l8| ) contribution. The accounted diagrams 
with dimension 9 operators have no loop integrations. There arc others v.e.v. of 
dimensions d < 9 particularly containing gluonic fields. All of them, however, 
correspond to at least one loop integration and are suppressed by the numerical 
factor (2tt)~ 2 . For this reason they are disregarded. 
The sum rule for X is given by 



C M 2 = -1 + ^e m2 / M2 \ a 2 L 4 / 9 + 6tt 2 f 2 M 4 E 1 ^ L^ 9 



9A 2 




U 1 M 2 I 4 M 2 9 Jo M 2 I 



(i 



Here M 2 is the Borel parameter, Ajy is defined as A^ = 327r 4 A 2 v = 2.1 GeV 
(0\r/\p) — \nv p , where v p is proton spinor, W 2 is the continuum threshold, W 2 = 
2.5 GeV 2 , 

a = -(27r) 2 (0|gg|0> = 0.55 GeV 3 (21) 

E (x) =l-e- x , E t (x) = 1 - (1 + x)e~ x 

L = ln(M / A) / ln(fi / A) , A = Aqcd = 200 MeV and the normalization point /i 
was chosen /i = 1 GeV. When deriving ( ^p| ) the sum rule for the nucleon mass was 
exploited what results in appearance of the first term, -1, in the right hand side (rhs) 
of (^0|). This term absorbs the contributions of the bare loop, gluonic condensate 
as well as a s corrections to them and essential part of terms, proportional to a 2 
and ttiqCl 2 . The values of the parameters, a, X%,W 2 taken above were chosen by 
the best fit of the sum rules for the nucleon mass (see [Q, Appendix B) performed 
at A = 200 MeV. It can be shown, using the value of the ratio 2m s /(m tI + 
m d ) = 24.4 ± 1.5 @ that o(l GeV) = 0.55 GeV 3 corresponds to m s (l GeV) = 
153 MeV. a s corrections are accounted in the leading order (LO) what results in 
appearance of anomalous dimensions. Therefore A has the meaning of effective A 
in LO. The unknown constant Co in the left-hand side (lhs) of (p^) corresponds to 
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the contribution of inelastic transitions p — ► -/V* — * interaction withA^ — > p (and 
in inverse order). It cannot be determined theoretically and may be found from M 2 



dependence of the rhs of (20) (for details see [|18|, |20J] ) . The necessary condition of 
the validity of the sum rule is |E| 3> \CoM 2 \exp[(— W 2 +m 2 )/M 2 ] at characteristic 
values of M 2 pc|]. The contribution of the last term in the rhs of (p0|)) is negligible. 
The sum rule (|20]) as well as the sum rule for the nucleon mass is reliable in the 
interval of the Borel parameter M 2 where the last term of OPE is small, less than 
10 — 15% of the total and the contribution of continuum does not exceed 40 — 50%. 
This fixes the interval 0.85 < M 2 < 1.4 Gey 2 .The M 2 -dependence of the rhs of 
( p0| ) at /q = 3 x 10~ 2 GeV 2 is plotted in fig. la. The complicated expression in rhs 
of (EC?) is indeed an almost linear function of M 2 in the given interval! This fact 
strongly supports the reliability of the approach. The best values of £ = E^ 1 * and 
Cq = Cq 11 are found from the y 2 fitting procedure 

n 

X 2 = - E - C tM? - R(M 2 )} 2 = mm, (22) 

i=l 

where R( 2 ) is the rhs of (p0|). 




Fig. 1. a) The M 2 -dependence of S + C M 2 at / 2 = 3 x 10" 2 GeV 2 ; b) E (solid 
line, left ordinate axis) and y%) (dashed line, right ordinate axis), as a functions 
of/o 2 - 



The values of E as a function of /q are plotted in fig. lb together with y X 2 - 
In our approach the gluonic contribution cannot be separated and is included in 



E. The experimental value of E can be estimated [£l],|22| (for discussion see 23 ) 
as E = 0.3 ± 0.1. Then from fig.lb we have / 2 = (2.8 ± 0.7) x 10~ 2 GeV 2 and 
X'(0) = (2.3 ± 0.6) x 10~ 3 GeV 2 . The error in /q and x' besides the experimental! 
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error includes the uncertainty in the sum rule estimated as equal to the contribution 
of the last term in OPE (two last terms in Eq.(pp|) and a possible role of NLO a s 
corrections. Allowing the deviation of \f~^ 1 by a factor 1.5 from the minimum 
we get x' = (24 ± 0.6) x 10~ 3 GeV 2 and E = 0.32 ± 0.17 from the requirement 
of selfconsistency of the sum rule. At /q = 2.8 x 10~ 2 GeV 2 the value of the 
constant Co found from the fit is Co = 0.19 GeV~ 2 . Therefore, the mentioned 
above necessary condition of the sum rule validity is well satisfied. 

Let us discuss the role of various terms of OPE in the sum rules ( |20] ) . To analyze 
it the sum rule ( |20| ) was considered for 4 different cases, i.e. when it is taken into 
consideration: a) only contribution of the operators up to d=3 (the term -1 and 
the term, proportional to f§ in (|2C|)); b) contribution of the operators up to d=5 
(the term ~ ho is added); c) contribution of the operators up to d=7 (three first 
terms in (po|)), d) the result (po|), i.e. all operators up to d=9. For this analysis the 
value of fo — 0.03 GeV 2 was chosen, but the conclusion appears to be the same 
for all more or less reasonable choice of f§. Results of the fit of the sum rules are 
shown in Table 1 for all four cases. The fit is done in the region of Borel masses 
0.9 < M 2 < 1.3 GeV 2 . In the first column the values of E are shown , in the second 
- values of the parameter C, and in the third - the ratio 7 = |-y/jp/S|, which is 
the real parameter, describing reliability of the fit. From the table one can see, 
that reliability of the fit monotonously improves with increasing of the number of 
accounted terms of OPE and is quite satisfactory in the case d 



Table 1. 



case 


E 


C{GeV- 2 ) 


7 


a) 


-0.019 


0.31 


lO" 1 


b) 


0.031 


0.3 


5.10~ 2 


c) 


0.54 


0.094 


9.10- a 


d) 


0.36 


0.21 


1.3 ■ 10~ 3 



3. Low energy theorems. 

Consider QCD with Nt light quarks, m ; <§; M ~ 1 GeV, i = 1, ...Nf. Define the 
singlet (in flavour) axial current by 

N f 

U$i x ) =^2ui( x h^7sq(x) (23) 

i 

and the polarization operator 

P^u (q)=i f d^xe^ (0 | T {^ 5 (x) , j, 5 (0) } | 0) . (24) 



The general form of the polarization operator is: 

P^(q) = -P L {q 2 )S^ + P T {q 2 )(S^q 2 + q^q v ) (25) 
Because of anomaly the singlet axial current is nonconserving: 

drf^ix) = 2N f Q 5 (x) + D(x), (26) 
where Qb{x) is given by (1) and 

N f 

D(x) = 2i^m i g i (x)7 5 (ji(x) (27) 

i 

It is well known, that even if some light quarks are massless, the corresponding 
Goldstone bosons, arising from spontaneous violation of chiral symmetry do not 
contribute to singlet axial channel (it is the solution of U(l) problem), i.e. to 
polarization operator P^ v {q). Pi(g 2 ) also have no kinematical singularities at q 2 = 
0. Therefore 

P„u{a)%qu = -P L (q 2 )q 2 (28) 

vanishes in the limit q 2 — > 0. Calculate the left-hand side (lhs) of ( p8| ) in the 
standard way - put q^qv inside the integral in ( pi| ) and integrate by parts. (For 
this it is convenient to represent the polarization operator in the coordinate space 
as a function of two coordinates x and y.) Going to the limit q 2 — > we have 



lim P^(q)W» = * / d x (° I T{2N f Q 5 (x), 2N f Q 5 (0) 

g 2 -»0 J 

+2N f Q 5 (x), D(0)+D(x), 2N f Q 5 (0)+D(x), D(0)} | 0) 
+4^mi(|? i (0) % (0) |0> 

i 

d 4 x(0 | [j 05 (x), 2N f Q 5 (0) } | 0)5(a:o) = (29) 

In the calculation of ( p9] ) the anomaly condition ( p6| ) was used. The terms, propor- 
tional to quark condensates arise from equal time commutator [705 (x), D(0)] XQ=0 , 
calculated by standard commutation relations. Relation ( |29[ ) up to the last term 
was first obtained by Crewther The last term, equal to zero according to stan- 
dard commutation relations and omitted in is keeped. The reason is, that we 
deal with very subtle situation, related to anomaly, where nonstandard Schwinger 
terms in commutation relations may appear. (It can be shown, that, in general the 
only Schwinger term in this problem is given by the last term in the lhs of (^9|) : no 
others can arise.) Consider also the correlator: 
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P M (q) =il d 4 x e« x {0 | T{j^ 5 (x), Q 5 (0)} | 0) (30) 

and the product P fJ _(q)q fl in the limit q 2 — > (or q 2 of order of the to 2 , where 
to w is the mass of Goldstone boson). The general form of P^q) is P^q) = Aq^. 
Therefore nonvanishing values of P^q^ in the limit q 2 — > (or of order of quark mass 
to, if | g 2 |~ to 2 - this limit will be also intersting for us later) can arise only from 
Goldstone bosons intermediate states in ([50]). Let us us estimate the corresponding 
matrix elements 

<0 | J M 5 | tt) = F Qll (31) 

(0 I Q 5 I tt) = F' (32) 

F is of order of m, since in the limit of massless quarks Goldstone bosons are 
coupled only to nonsinglet axial current. F' is of order of m 2 /^ ~ to, where / w is 
the pion decay constant (not considered to be small), since in massless quark limit, 
the Goldstone boson is decoupled from Q5 . These estimations give 

q 2 

P/j,q» — 2 1™ 2 ( 33 ) 

q A — m£ 

and it is zero at q 2 — > 0, to 2 ^ and of order of to 2 at q 2 ~ to 2 . In what follows I 
will restrict myself by the terms linear in quark masses. So, I can put P ll {q)q ll = 
at q — > 0. The integration by parts, in the right-hand side (rhs) of (|30|) gives: 



lim = - / cf x(0 | T{2N f Q 5 (x), Q 5 (0) + D(a:), Q 5 (0)} | 0) 



d 4 x(0 I [ j 05 (x), Q 5 (0) ] I 0)5(x ) = (34) 
After the substitution of |J]3n ( p9[ ) arise the low energy theorem: 

1 f d 4 z(0 I T{2iV / Q 5 (x), 2JV>Q 6 (0)} | 0) 



N f 

d 4 x(0 I T{D(x), D(0)} I 0) -4^to,(0 | ?i(0) ?i (0) | 0) 

i 

+ 1 f d 4 x(0 I [ j$ 5 (x), 2N f Q 5 (0) } I Q)8(xo) = (35) 



The low energy theorem (|35|), with the last term in the lhs omitted, was found by 
Crewther [||. 
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4. One and two light quarks. 

Consider first the case of one massless quark, Nf — 1, m = 0. This case can easily 
be treated by introduction of #-term in the Lagrangian, 

AL = 8^- G%G^ V (36) 

The matrix element (0 | Q5 | n) between any hadronic state | n) and vacuum is 
proportional 

d 4 x(0 I Q 5 (x) I n) ~ (0 I ^InZ \ n)e=o, (37) 

where Z = e lL and L is the Lagrangian. The gauge transformation of the quark 
field ip' — ► e M7£ "0 results to appearance of the term 

SL = adfj^ = a(a s /47r)G^G^ (38) 

in the Lagrangian. By the choice a = —6 the 0-term ( ^6| ) will be killed and 
(d/d6)lnZ = 0. Therefore, x(0) = (Crewther theorem). The first term in @ 
vanishes, as well as the second and third, since m = 0. From ( |35| ) we have, that 
indeed the anomalous commutator vanishes 

(0| [jas(x), Q 5 (0) ]x =o I 0) = 0, (39) 

supporting the assumptions done in 

Let us turn now to the case of two light quarks, u,d,Nf = 2. This is the case 
of real QCD, where the strange quark is considered as a heavy. Define the isovector 
axial current 

j% } = (ui„7 5 u - d W d)/V2 (40) 
and its matrix element between the states of pion and vacuum 

(0 I S I = ( 41 ) 

where q M is pion 4- momentum, f n = 133 MeV. Multiply ( (fl| ) by q^. Using Dirac 
equations for quark fields, we have 

2i - i 

~y=(° I m u uj 5 u - m d dj 5 d \ ir) = — (0 | {m u + m d )(u-f 5 u - dj 5 d) 

+ (m« - m d )((uj 5 u + d'ysd) I 7r) = /^mj, (42) 

where m u , m d are u and <i quark masses. The ratio of the matrix elements in lhs 
of (143) is of order 
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(0 | wy 5 u + d<y 5 d | tt) ^ m M - m d 

(0 | u 75 u - J 76 d | tt) ~ M ' 1 j 

since the matrix element in the numerator violates isospin and this violaton (in 
the absence of elecntomagnetism, which is assumed) can arise from the difference 
m u — md only. Neglecting this matrix element we have from (ptj 



4=<0 | u l5 u - d l5 d | tt) = hr f* (44) 
V2 m u + m d 



Let us find x(0) from low energy sum rule ( p5[ ) restricting ourself to the terms linear 
in quark masses. Since D(x) ~ m, the only intermediate state contributing to the 
matrix element 

J d A x(i) | T{D(x), D(0)} | 0) (45) 
in (35) is the one-pion state. Define 

D q = 2i(m u u"f 5 u + m d d~i 5 d). (46) 



Then 



(0 | D q | tt) = i(0 | (m u + m d ){u"f 5 u + dj 5 d) + (m u - m d ){uj 5 u - d-y 5 d) \ tt) 



= V2 Uml, 47 

m u + m d 

where the matrix element of singlet axial current was neglected and ( ^i|) was used. 
The substitution of (f47| ) into (ji^) gives 



i f d^e^iO | T{D q (x), D q (0)} | 0)^ = lim | -^J-_ 2 f m " md ) f*m* 
J q^o I q* -m* ym u + m d J 

= _ 4 (mu-m d r 
m u + m d 

In the last equality in (Q) Gell-Mann-Oakes-Renner relation 

1 f 2 m 2 

(0 I qq 10) = -- \" (49) 

2 m u + m d 

was substituted as well the SU (2) equalities 

(Q\uu\0) = {0\dd\0} = {a\qq\ 0). (50) 
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From @ and (g|) we finally get: 

X (0)=i [ d*x(0\T{Q 5 (x), Q B (0)}|0) 



in concidence with eq.3. 

In a similar way matrix element (0 | Q5 \ -k) can be found. Consider 



(o I m I o> (5i) 



The estimation of F gives 



(0 I J M 5 I 7T) = Fg„ 



~ M ^ 



(52) 



(53) 



and after multiplying of (|32|) by <? M the rhs of ( p2| ) can be neglected. In the lhs we 
have 



The substitution of 



(0 I D q I tt) + 2JV>(0 I Q 5 I tt) = 
into (|5l) results in 



(0 I Q 5 



1 IUi 



(54) 



(55) 



2^2 w„ + m d 

The relation of this type (with a wrong numerical coefficient) was found in [pl|, the 
correct formula was presented in p5|] . From comparison of ( |5~i| ) and (|55|) it is clear, 
that it would be wrong to calculate x(0) by accounting only pions as intermediate 
states in the lhs of (|5l]) - the constant terms, reflecting the necessity of subtraction 
terms in dispersion relation and represented by proportional to quark condensate 
terms in ( |35"| ) are extremingly important. The cancellation of Goldstone bosons pole 
terms and these constant terms results in the Crewther theorem - the vanishing of 
x(0), when one of the quark masses, e.g. m u is going to zero. 



5. Three light quarks. 

Let us dwell on the real QCD case of three light quarks, u, d and s. Since the 
ratios m u /m s , md/m s are small, less than 1/20, account them only in the leading 
order. When the u and d quark masses m u and md are not assumed to be equal, 
the quasi-Goldstone states w° and 77 are no more states of pure isospin 1 and 
correspondingly: in both of these states persist admixture of other isospin ]lO|,pj}| 
proporional to m u — md- (The violation of isospin by electromagnetic interaction 
is small in the problem under investigation [jl0| and can be neglected. The 77' — r\ 
mixing is also neglected.) In order to treat the problem it is convenient to introduce 
pure isospin 1 and pseudoscalar meson fields t/?3 and ip% in SU (3) octet and the 
corresponding states | P 3 ), | P 8 ) Then in the SU(3) limit 



(0 I j% I Ps) - Uq„ 



(56) 
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(0 I J$ I Ps) = Uq» (57) 

where 



3 J = ("7 M 75W + ^75^ - 2s7 Al75 s)/%/6 (58) 



and j'^g is given by (|40|). The states | P3), | Pg) are not eigenstates of the Hamil- 
tonian. In the free Hamiltonian 

H = -mlifj + -m\<p\ + (P 8 I P 3 )</W8 + kinetic terms (59) 

the nondiagonal term ~ (f3(ps is present ( to \ and ml in ( p9| ) coincides with to^ 
and to^ up to terms quadratic in (Pg | P3)). The nondiagonal term was calculated 
in (!(]] on the basis of PC AC and current algebra (see also [p] ) 

, „ , „ , 1 9 TO,, — TOrf 

P 8 P 3 = -=to" " 60 

V3 TO„ + TO d 

The physical ir and 77 states arise after orthogonalization of the Hamiltonian ( |59| ) 

I tt> = cosO I P 3 ) - sin6 I P 8 ) 

I 77) = sine | P 3 ) + cos6» | P 8 ) (61) 
where the mixing angle e is given by (at small 6) 



e _ (Ps I P3) ^ (Ps I Pi) = J_ m| m M -TO d 

TO^ - TO.2 TO^ ^ TO.2 TO U + m d 

In terms of the fields ^3 and ^8 PCAC relations take the form llOi 



(62) 



W$ = U«V3 + (Ps I PM) (63) 

dJS = + ( p s I Pi)<Pz) (64) 

Our goal now is to calculate the contribution of pseudoscalar octet states to the 
second term in the lhs of (|35|). It is convenient to use the full set of orthogonal 
states I P3), I Pg) as the basis. Use the notation 

D = D q + D s D s = 2im s s"f 5 s, (65) 
where D q is given by (p3|). The matrix element 

(0 I D g I P 3 ) = V2 mu ~ md Uml (66) 
m u + m d 

can be found by the same argumentation, as that was used in the derivation of ( |4^ ) . 
In order to find (0 | D q | Pg) take the matrix element of eq.|6^ between vacuum and 
\Pr) 
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(67) 



The substitution in the lhs of (|67j) of the expression for c^iig through quark fields 
gives 



(0 | uj 5 u - dj 5 d | P$) 



m u - m d 
m u + m d 



(0 | u lb u + d lb d | Pft)- 



. 2 2 m u -m d 
V 3 (m u + m d y 



(68) 



In a similar way take matrix element of eq. (|64|) between (0 | and | Ps) and 
substitute into it (p8|). We get 



V6 



(0 



m u + m d 



(m u - m d ) 2 
m u + m d 



{u-f b u + dj 5 d) - 4m s S7 5 s | P s ) 



,2 2 1 2 , ( m u ~ m d 

3 \m u + m d 



(69) 



As follows from SU(3) symmetry of strong interaction 

(0 | u lb u + d l5 d | P s ) = -(0 | s l5 s | P 8 ) (70) 
up to terms of order m q /M, which are neglected. From (|69|),(|7C|) we find: 



(0 | D s | P a ) = - 



_ 1 (m u - m d ) 2 
4 m s (m u + m d ) 



(Q\D q \P 8 )=4d-Umi 



1 



m u m d 



m s (m u + m d ) 



m u m d 



V 3 (m u + m d y 
in notation (B^),(^). When deriving (|7l|) the .St/ (3) relation 



2 4 2 
m r, = o m 7r 



(71) 



(72) 



(73) 



was used. In (71) the small terms ~ m u /m S: m d /m s are accounted, because they 
are multiplyed by large factor mi. In (f72|) small terms are disregarded. The matrix 
element (0 | D s \ P 3 ) can be found from (|64|). We have 



-j=(0\D g -2D s \P 3 )=U(P s \P 3 ) 



(74) 



The substitution of <M) and (BOh into (WW) gives 
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(0 | D s | P 3 ) = (75) 

Equations (|66|) , ([n]) , (|72]) and ((75]) allow one to calculate the interesting for us cor- 
relator 

i J d 4 a;(0 | T{D(x), D(0)} | 0) (76) 

when the sets | P3KP3 | and | Ps)(Ps | are taken as intermediate states. But 
I P3), I Pg) are not the eigenstates of the Hamiltonian, they mix in accord with 
([59l). Therefore the transitions (P$ | P3) arising from the mixing term in (^9|) must 
be also accounted. There are two such terms. The one corresponds to the transition 
(0 I D s I P 8 ) (P 3 I D q I 0) and its contribution to @ is given by 

liin (-2(0 I D s I P 8 )—±—( Ps I p 3 )^— { p 3 \ Dq \ ) 



*o [ ' q 2 — mi q 2 — 

= V>l( r ^-^) 2 (77) 
y m u + m d J 

The other corresponds to the transition between two D s operators, where (P3 | P3) 
enter as intermediate state. This contribution is equal to: 



lim {-(0\D S I P 8 ) 2 1 2 (P 8 I P 3 ) 2 1 2 (P 8 I P 3 ) 2 1 2 (P 8 I D s I 0) 




(78) 

It is enough to account only matrix elements, with D s operators, since they are 
enhanced by large factor m 2 . All others are small in the ratio m^/rn 2 . 
Collecting all together, we get: 



i / d 4 x(0 I T{D{x), D(0)} I 0) = flml { 2 



m u - m d 



m u m d 



m u + m d I (m„ + m d ) 2 



+2- 



1 m u + m d 



m u + m d 



+2 



1 \ 2 / v 2 

m u - m d \ 1 m u ~ m d \ 



1 m u - m d 

2 m s (m u + m d ) 

2m s 



1 - 2- 



m u m d 



m s (m u + m d ) 



j-2 2 



m u + m d 2 \ m u + m d 



9 / m u - m d \ 5 

m u +m d 2\m u +m d J 2 

(79) 
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The first term in the figure bracket in ( |79| ) comes from (0 | D q | P3) 2 , the second - 
from (0 I DA P 8 ) x (P 8 \ D s | 0V, the third - from (0 | D s \ P 8 ) 2 , the last two terms 
are from (|77|) , (wm . Adding to (f79|) the proportional to quark condensate term 

4(m u +m d + m s )(0 \ qq | 0) (80) 



in (35), we finally get for 3 quarks at m Ul md <C m s 



1 / d^O I T{2N f Q 5 (x), 2N f Q 5 (x)} | 0) = 36 m " m " (0 | qq | 0) (81) 
' m u + m d 



and 



X(0)= ^ <0 I gg I 0>, (82) 

since in this case ANj = 36. Eq.82 coincides with eqj|[ obtained in N c — * co 
limit. This fact demonstrates, that iV c — > 00 limit is irrelevant for determination 
of x(0) (at least for the cases of two or three light quarks and at m u ^m d <§; m s ). 
x(0) for three light quarks at m u , m d <C m s - eq.p2 coincides with x(0) in the two 



light quark case, [see |8j] and (51)] i.e in this problem, when m u ,md <§C m s , there is 
no difference if s-quark is considered as a heavy or light - it softly appears in the 
theory. 

Determine the matrix elements (0 | Q5 \ ?/) and (0 | Q5 | 7r). Following ]26| 
consider 

(0 I J M 5 I V) = (83) 

F is of order of f 7T (m s /M) and can be put to zero in our approxamation. By taking 
the divergence from (^), we have 

(0 I D s + 6Q 5 I 77) = (84) 

The use of ( |si| ) gives (the ir — r\ mixing as well terms of order m u /m s , md/m s may 
be neglected here): 

(0 I Q 5 I V) = \\j\Um 2 n (85) 
Relation ( ^5| ) was found in [^6| . By the same reasoning it is easy to prove that 

(0 I D q I P 3 ) + (0 I D s I P 3 ) + 6(0 I g 5 I P 3 ) = (86) 



The first term in (|36|) is given by (|66| ), the second one is zero according (|75j). For 
the last term we can write using 

<0 I Q 5 I P 3 ) = <0 I Q 5 I tt) + 9(0 I Q 5 I P 8 > (87) 
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Eq.'s (^,(|87|) give 

, , „ , , 1 o m„ — rrLri 

<PI*l->— < 88 > 

- the same formula as in the case of two light quarks. 

It is clear, that the presented above considerations can be generalized to the 
case, when u,d and s-quark masses are comparable. The calculation became more 
cumbersome, but nothing principially new arises on this case. 



6. q 2 — dependence of x(g 2 ) at low q 2 . 

Let us dwell on the calculation of the q 2 -dependence of x(l 2 ) at l° w I Q 2 I m QCD, 
restricting ourselves by the first order terms in the ratio q 2 /M 2 , where M is the 
characteristic hadronic scale, M 2 ~ 1 GeV 2 . In this domain of q 2 xil 2 ) can be 
represented as 

X(q 2 ) = X(0) + X'(0)q 2 + R(q 2 ) - R(0) (89) 

x(0) for the QCD case-three light quarks with m u ,m,d <S m s - was determined 
in Sec.V. x'(0) (its nonperturbative part) for massless quarks was found in |lj 
basing on connection of x'(0) with the part of proton spin carried by u,d,s 
quarks. Its numerical value is given by ([!]). What is left, is the contribution of light 
pseudoscalar quasi-Goldstone bosons R{q 2 ), which has nontrivial g 2 -dependence 
and must be accounted separately. R{q 2 ) vanishes for massless quarks and did not 
contribute to x'(Q), calculated in H. i?(0) must be subtracted from R(q 2 ) since it 
was already accounted in x(0). R(q 2 ) can be written as 

R(q 2 ) = -(0 | Q 5 I t) 2 ^^ - (0 | Qb I V) 2 ^-^ (90) 
q* — q z — 

The problem of rj — ir mixing is irrelevant in the difference R(q 2 ) — R(0) in any 
domain Ig 2 |~ m 2 and | q 2 |~ m 2 . The matrix elements entering (|9(]) are given 
by (|S5l) , (|8^) . Taking the difference R(q 2 ) — R{0) and using the Euclidean variable 
Q 2 = —q 2 , we have 



x(Q 2 )=x(0)-x'(0)Q 2 --f 2 Q 2 



m d 



3Q 2 



(91) 



Eq.(^Tj) is our final result, where x(0) is given by © and X '(0) by §). The 
accuracy of (|9l| ) is given by the parameters Q 2 /M 2 , m 2 /M 2 , m 2 /M 2 <§C 1, (two last 



characterize the accuracy of SU(3) x SU(3)). At Q 



the last term comprise 



about 20% of the second (the first term is very small, x(0) ~ —4 ■ 10~ 5 GeV 4 ). 
Evidently, the last term is much bigger in the Minkovski domain, Q 2 < 0, since 
there are pion and r\ poles. As was mentioned above, x(0) found here concides with 
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X(0) obtained in by considering large N c limit. However, the <? 2 -dependence is 
completely different. Namely, for x(<7 2 ) m was found the relation (eq.(A4') in 
@) 

(92) 

where the Goldstone boson masses /if are related to quark condensate by 

/i 2 = ~2m~(0 \qq\0) (93) 

and a is some constant of order of hadronic mass square. At q 2 = follows cq. pj fo r 
x(0) if the inequality a/N c fj, 2 » 1 is assumed. However, at \q 2 \ ~ m 2 , m 2 (|92|) 
strongly differs from (|9l|): in ( p^ ) there are zeros at the points q 2 = ml , m 2 , but 
not poles, as it should be and as it take place i n (pl| ). And also the most important 
at low Q 2 hadronic term x'(0)Q 2 is absent in (|92[). 



X(l 2 ) = 



aFl 
2N r 



1 - 



N c ^q 2 - 



7. Summary. 

The (^-dependence of topological charge density correlator x(<Z 2 ) (2) in QCD was 
considered in the domain of low q 2 . 

The QCD sum rule, connecting x'(0) (for massless u, d, s quarks) with the part 
of the proton spin S, carried by u,d,s quarks in the polarized proton, is derived. 
The numerical value of x'(0) was found from the requirement of selfconsistancy of 
the sum rule. In the limit of errors the same value of x'(0) follows also from the 
experimental data on E. 

For the cases of two and three light quarks the values of x(0), obtained earlier 
0"|§ were rederived basing on the low energy theorems and accounting of quasi- 
Goldstone boson (jr, rj) contributions. No large N c limit was used and it was no 
appeal to the 0-dependence of QCD Lagrangian (except of the proof of absence of 
anomalous commutator in the sum rule (|35|)). The only concept, which was used, 
was the absence of Goldstone boson contribution as intermediate state in the singlet 
(in flavour) axial current correlator in the limit of massless quarks. In the three 
light quark case - the case of real QCD - the mixing of 7r and rj is of importance 
and was widely exploited. The q 2 dependence of x(9 2 ) was found as arising from 
two sources: 

• 1. The contribution of hadronic states (besides 7r and 77), determined from 
the connection of x'(0) with the part of the proton spin, carried by quarks. 

• 2. The contributions of 7r and 77 intermediate states. These contributions were 
calculated by using low energy theorems only. The final result is presented in 



eq.91 
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